Volume 94, Number 2, March-April 1989 

Journal of Research of the National Institute of Standards and Technology 



A Supercritical Fluid Chromatograph 
for Physicochemical Studies 



Volume 94 



Number 2 



March-April 1989 



Thomas J. Bruno 

National Institute of Standards 
and Technology, 
Boulder, CO 80303 



A supercritical fluid chromatograph has 
been designed and constructed to make 
physicochemical measurements, while 
retaining the capability to perform 
chemical analysis. The physicochemical 
measurements include diffusion coeffi- 
cients, capacity ratios, partition coeffi- 
cients, partial molar volumes, virial 
coefficients, solubilities, and molecular 
weight distributions of polymers. In this 
paper, the apparatus will be described in 
detail, with particular attention given to 
its unique features and capabilities. The 



instrument has recently been appUed to 
the measurement of diffusion coeffi- 
cients of toluene in supercritical carbon 
dioxide at a temperature of 313 K, and 
pressures from 133 to 304 bar 
(13.3-30.4 MPa). The data are dis- 
cussed and compared with previous 
measurements on similar systems. 
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1. Introduction 



The methods of chromatography have been ap- 
plied to physicochemical problems such as thermo- 
physical property determination for nearly 30 years 
[1,2]. This application of chromatography to non- 
analytical problems stems from an understanding of 
the physical and chemical processes which are 
known to occur during chromatographic separa- 
tions [3]. As an example, since hydrogen bonding 
can play a role in chromatographic separation, we 
may apply chromatography to the study of hydro- 
gen bonding thermodynamics [4]. 

The development of chromatography during the 
last 80 years can be divided into distinct historical 
periods [5], each with its innovations, fads and fail- 
ures. During the current period, we have seen the 
emergence of supercritical fluid chromatography 
(SFC, in which the carrier is a fluid held above its 
critical point). Some properties of a typical super- 
critical fluid can be seen in table 1. The density of 



Table 1. Comparison of representative fluid properties 



Gas Liquid Supercritical fluid 



Density, p 10"^ 

g/mL 

Diffusivity, D 10"' 

cmVs 

Dynamic viscosity, tj 10"* 
g/(cm-s) 



5X10- 



10- 



0.7 



10- 



10- 



the supercritical fluid is very similar to that of a 
liquid phase. This property explains the greatly en- 
hanced solvation power of the supercritical fluid 
with respect to the gas phase. The viscosity of the 
supercritical phase closely resembles that of the gas 
phase, thus allowing for easy mass transfer. The 
thermal conductivity (though not shown in table 1) 
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is also relatively large, as one would expect from 
density and viscosity considerations. The diffiisiv- 
ity (self-diffusion) of the supercritical phase is inter- 
mediate between that of a gas and a liquid. This 
property gives the supercritical fluid the advantage 
over liquid-liquid extraction. There are many ex- 
cellent reviews describing the advantages and ap- 
plications of SFC, and the reader is referred to 
these for additional details [6,7]. 

Most of the applications of SFC found in the 
literature involve analytical or separation prob- 
lems. The application to physicochemical studies 
has been relatively slow due to many experimental 
difficulties. Nonetheless, SFC has been applied to a 
number of thermophysical problems by several 
groups. This has included the measurement of ca- 
pacity ratios, partition coefficients, binary diffusion 
coefficients, partial molar volumes, virial coeffi- 
cients, solubilities, and polymer molecular weight 
distributions. In this work, the physicochemical su- 
percritical fluid chromatograph has been applied to 
the measurement of binary diffusion coefficients 
[8-16]. This diffusion coefficient describes the ten- 
dency of the solute to diffuse into the carrier (usu- 
ally referred to as the solvent). 



2. Theory 

Chromatographic instruments and methodology 
has been applied to the measurement of binary dif- 
fusion coefficients since the early 1960s [17,18]. 
The theoretical principles of the so-called peak 
broadening method were formalized by Taylor and 
Aris [19-23]. Since many authors have treated the 
theory in full detail, only a brief summary of the 
main results will be provided here. 

A sharply defined spike of solute (8-function 
like) introduced into a laminar stream of a carrier 
flowing in an uncoated tube will be subject to both 
convective action along the axis of the tube and 
molecular diffusion in the radial direction. If the 
tube is assumed to be straight, the mathematical 
treatment of Taylor and Aris gives a simple expres- 
sion for the plate height H: 



hJ^^ 



24D, 



(1) 



where Dn is the binary diffusion coefficient of the 
solute into the carrier (solvent), u is the linear ve- 
locity of the mobile phase solvent, and r is the in- 
ternal radius of the uncoated tube. The plate height 
jy is an important and experimentally accessible 



quantity in chromatography, since it describes the 
efficiency of the chromatographic system [24]. 
This quantity is the width of a peak (as designated 
by its variance, cr^, in length units as opposed to 
time units) relative to the distance traversed inside 
the column or tube (i.e., the length of the tube, L): 



H = cr'/L 



(2) 



In a straight tube, the concentration profile of the 
solute in the carrier will become Gaussian-like 
when jy<0.02 m [25]. The same result may be ob- 
tained from the Golay equation (mathematical de- 
scription of capillary gas chromatography) when 
one assumes a zero thickness for the stationary 
phase film (i.e., an uncoated tube) [26]. At the car- 
rier fluid velocities encountered in practice, eq (1) 
reduces to: 



H = 



24Dn ' 



(3) 



which allows calculation of the binary diffusion 
coefficient. 



3. Experimental 

The supercritical fluid chromatograph con- 
structed for this work is shown schematically in 
figure 1. The main components of the instrument 
are the solvent delivery system, the sample injec- 
tion system, the column thermostat, and the detec- 
tion/quantitation system. Safety devices to prevent 
column temperature and pressure runaway, and in- 
ert gas purging of key components provide for ex- 
plosion-proof operation. 

The solvent delivery system is designed to han- 
dle fluids which are either gaseous or liquid at 
room temperature and pressure. The pump is a 
modified commercial, double piston pressure-con- 
trolled device which is electronically compensated 
for pressure pulsation. This pump is capable of de- 
livering a pressure of 100 MPa, and all associated 
transfer lines are rated to safely accommodate this 
limit. The pump head is enshrouded by a welded- 
seam stainless steel (304) cover (0.16-cm thickness), 
the inner surface of which is insulated with a 1.27- 
cm layer of foamed silica. Brazed to the front sur- 
face of the cover is a 0.79-cm (inside diameter) 
compression fitting which accepts the cold stream 
of a Ranque-Hilsch vortex tube. The vortex tube is 
used to cool the pump head and carrier stream to 
facilitate delivery of lower critical point fluids such 
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Figure 1. Schematic diagram showing the major components of the supercritical fluid chromatograph used in this work. 
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as carbon dioxide [27,28]. Before the carrier fluid 
enters the pump head, it is first passed through a 
heat exchanger which consists of a 600-cm length 
of stainless steel (304) tubing (0.076-cm inside di- 
ameter, 0.16-cm outside diameter). This heat ex- 
changer insures that fluid is delivered to the pump 
head as a liquid. The vortex tube is operated in an 
intermediate mode (between maximum cooling and 
maximum temperature difference, with an applied 
air pressure of 0.7 MPa), and provides a carrier 
stream temperature of —20 °C. This minimizes va- 
por locking and cavitation inside the pump heads. 
The pump is followed by a pulse dampener and 
pressure transducer. The pulse dampener is a coil 
of flattened stainless steel tubing (0.64-cm outside 
diameter) which absorbs much of the low level pul- 
sation not handled by the electronic compensation. 
This component is necessary since no piston pump 
can operate in true pulse-free fashion. The pressure 
transducer is a strain gauge device which is cali- 
brated periodically using a high precision Bourdon 
tube transfer standard. This transfer standard is it- 
self calibrated using a dead weight pressure balance 
traceable to the NIST primary standard. The un- 
certainty in the measured pressure has been deter- 
mined at ±0.40 bar (±0.040 MPa). 



After leaving the pressure transducer the fluid 
enters a heat exchanger (fig. 2) inside the column 
oven, where the flash to supercritical temperature 
occurs. This heat exchanger is a 300-cm section of 
stainless steel (304) tubing (0.32-cm outside diame- 
ter, 0.07-cm inside diameter). A vibrating tube den- 
simeter is downstream from the heat exchanger, to 
allow independent density measurements of the 
carrier if desired. The densimeter places operating 
constraints upon the entire system (41 MPa maxi- 
mum pressure, 160 °C maximum temperature), and 
must be removed for higher temperature or pres- 
sure operation. 

Following the solvent delivery system is the 
sample injector. The injector used in this work is 
the flow-through extractor coupled with a high- 
pressure chromatographic sampling valve shown 
in figure 3 [29]. This arrangement is most satisfac- 
tory for experiments involving repetitive injections 
of the same solute. An aliquot of solvent-borne 
sample is syringe-deposited into the extractor. The 
extractor is then heated and evacuated to remove 
the sample solvent. After the solvent is removed 
(an operation which takes approximately 5 min- 
utes), the extractor is filled with the supercritical 
carrier (solvent) at the same temperature and pres- 
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Figure 2. Diagram of ttie ttiermostatted column area containing the diffusion tube and vibrating tube densimeter. 
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Figure 3. Schematic diagram of the sample injection system. 



sure as that in the column. The sample is then 
loaded into the sample loop to allow injection into 
the carrier stream. The injection is achieved auto- 
matically under computer control using an pneu- 
matic actuator (with helium as the working fluid) 
equipped with pilot valves and a large ballast vol- 
ume to provide fast switching. 

Upon injection, the carrier-borne solute is trans- 
ported to the column area (fig. 2) inside of a modi- 
fied commercial forced-air oven capable of 
maintaining a temperature of 350 °C. The major 
modifications to the oven include baffles which 
promote uniform air movement and an inert gas 
purge line for safety. The column area contains the 
heart of the physicochemical experiment, which 
may consist of a coated or uncoated capillary or a 
packed column. In the case of diffusion coefficient 
measurements, the column consists of a long un- 
coated capillary. This capillary is currently a 3040- 
cm long continuous section of 316 stainless steel 
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tubing (0.159-cm outside diameter, 0.025+0.0013- 
cm inside diameter) which is connected directly to 
the injection valve. The length of the diffusion tube 
was determined by weight, using a calibration 
equation. This tube is coiled in a 30-cm diameter, 
held inside of an aluminum racetrack which inte- 
grates out local temperature variations. The race- 
track also serves as the support for the vibrating 
tube densimeter referred to earlier. Temperature 
measurement is provided by a platinum resistance 
probe located in the center of the racetrack. Six 
pairs of gradient-sensing thermocouples (type-j, 
thermally tempered) referenced to the main ther- 
mometer provide an indication of temperature 
nonuniformity. This nonuniformity is then reduced 
to a negligible level using a set of manually con- 
trolled low power shimming heaters. The race- 
track is supported inside the oven by titanium rods 
of 0.64-cm diameter. The low thermal conductivity 
of the titanium limits heat transfer into or out of the 
oven. The temperature uniformity of the racetrack 
has been measured at ±0.015 °C at temperatures at 
or above 40 °C. 

After leaving the racetrack, the solute is trans- 
ported to the detector (fig. 4) by the diffusion tube, 
as shown in figure 2. The detector is a modified 
chromatographic flame ionization detector con- 
tained in a separate oven directly above the column 
oven. The major modifications made to the com- 
mercial unit were in the sample and jet gas inlet 
lines. Provisions have also been made for the intro- 
duction of a make-up gas where needed. A fused 
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Figure 4. Schematic diagram of the modified flame ionization 
detector. 



quartz capillary restrictor, attached to the diffusion 
tube, releases the solute and carrier directly into 
the flame. The detector is operated to provide an 
approximate sensitivity of 10~" mol/s. The inside 
diameter and length of this restrictor capillary is 
chosen so as to maintain the carrier velocity be- 
tween 2 and 6 cm/s. The temperature of the detec- 
tor is maintained at 300 °C to prevent the carrier 
from cooling (and possibly solidifying) upon de- 
compression. This cooling action has been noted as 
a cause of baseline "spiking" in unheated detectors. 
The output from the detector is logged on a com- 
mercial electronic integrator, from which the re- 
tention times and peak widths may be extracted. 



4. Results and Discussion 

The main sources of error in this experiment 
stem from the diffusion tube itself, the sample injec- 
tion process, pressure drop across the tube, and ad- 
sorption of the solute on the tube walls. A brief 
discussion of these errors will be represented here, 
especially with respect to this apparatus. A number 
of excellent general reviews are available in the lit- 
erature [30-34]. 

The theory of Taylor and Aris, summarized ear- 
lier, is strictly based on straight tubes of circular 
cross-section and uniform inside diameter. For 
practical reasons the tube is held inside the oven as 
a helical coil rather than a straight length of tubing. 
In addition, tubes of perfectly circular cross-sec- 
tion and uniform inside diameter are an idealization 
not available in the laboratory. We must therefore 
determine to what extent our experimental appara- 
tus departs from the theoretically assumed condi- 
tions. It can be shown [34] that the coiling of the 
diffusion tube will have no harmful effects if cer- 
tain criteria are met. The ratio (o)) of the radius of 
the coil to that of the inside diameter of the diffu- 
sion tube should be greater than 100. In the present 
apparatus, the radius ratio is approximately 1200. 
Another requirement is that the inequality 
De^Sc <20 must be satisfied, where De is the Dean 
number and Sc is the Schmidt number (see Ap- 
pendix 1 for definitions). In the design of the 
present apparatus, De^Sc is always less than 10. 
Secondary flow effects are minimized by maintain- 
ing the carrier velocity below 6 cm/s, thus assuring 
flow in the laminar regime. 

The problem on nonuniformity of the tubing in- 
ternal radius is somewhat more difficult to control. 
It is possible to make corrections for this nonuni- 
formity if the nature of the departure is known 
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(for example, a sinusoidal variation superimposed 
on the radius has been considered [34]). Selection 
of a good quality seamless tubing will help mini- 
mize most of the problems associated with tube 
nonuniformity. It should be noted that the uncer- 
tainty in the tube's inside diameter cited in the 
experimental section (~5 percent) is a worst case 
limit, but this value has been used in the overall 
error analysis. Since the tube length and radius are 
both temperature dependent, corrections are ap- 
plied to account for the effect of thermal expansion 
on these parameters. The effect of the applied pres- 
sure on the internal tube radius is negligible. The 
effects of connections between the tubes have been 
minimized by design (only one such connection is 
used, having the same internal radius as the diffu- 
sion tube), and are considered to be negligible. 

The errors associated with injection have been 
addressed in several ways. The sample is intro- 
duced at the same pressure and temperature as that 
of the carrier stream, and at infinite dilution. The 
volume of the sample loop is very small (s;3.0 juL) 
as compared with the volume of the diffusion tube. 
The sampling loop is switched into the carrier 
stream only for a short time, and then switched 
back to the fill position. This has been found to 
decrease problems from solute adsorption. The in- 
jection process is done extremely fast using the pi- 
lot valve system, resulting in a negligible pressure 
pulse. 

The pressure drop across the diffusion tube has 
been measured at between 0.3 and 0.6 bar (0.03- 
0.06 MPa). Since this is on the order of the experi- 
mental error of the pressure measurement, no 
corrections are made, but close proximity of the 
critical point is avoided. Measurements made 
within a few degrees of the critical temperature 
and near the critical pressure will require consider- 
ation of this density gradient effect. Solute adsorp- 
tion on the inside walls of the tube has not been a 
problem in the present work, as judged from the 
peak symmetry. Adsorption is often a problem 
when using highly polar solutes at relatively low 
carrier density and temperature. To address this is- 
sue, diffusion tubes which have larger internal radii 
are used so as to decrease the surface area-to-vol- 
ume ratio. This must be done with consideration of 
the trade-off in radius ratio o), and the secondary 
flow effects which may result. 

As an example of the operation of this apparatus, 
measurements of binary diffusion coefficients of 
toluene (at infinite dilution) in supercritical carbon 
dioxide are presented in table 2 [35]. The measure- 
ments were made at 313.83+0.02 K, and at pres- 



sures from 133 to 304 bar (13.3 to 30.4 MPa). The 
carrier fluid density corresponding to these tem- 
perature-pressure pairs was calculated from the 
32-term Benedict- Webb-Rubin (BRW) equation of 
state for carbon dioxide [36], and ranged from 
0.746 to 0.910 g/cm^. As discussed earlier, the raw 
chromatographic data obtained were the peak 
widths and breakthrough times (the term "reten- 
tion time" being considered inappropriate due to 
the absence of a stationary phase). At each density, 
15 separate determinations were made, furnishing 
the experimental uncertainty for the error analysis. 
The combined uncertainties of all measured quanti- 
ties provide an overall estimate of between 5 and 6 
percent for the data presented here. The diffusion 
coefficients are all on the order of 10"" cmVs, and 
decrease with increasing carrier density as can be 
seen in figure 5. These data fit in quite well with 
previous data on lower molecular weight aromat- 
ics, although most of the previous data were taken 
at somewhat lower densities [8-11,14]. Compari- 
sons of this data with several predictive approaches 
are as yet incomplete and will be presented in the 
future. Current work also includes a study of a 
homologous series of straight chain hydrocarbons, 
and several members of the carotene family. 



Table 2. Measured binary diffusion coefficients, Dn, of toluene 
in supercritical carbon dioxide, at 313.83±0.02 K 



CO2 density 
(g/mL) 



A2XIO* 
cmVs 



0.746 
0.801 
0.893 
0.867 
0.890 
0.910 



1.30 
1.23 
1.21 
1.20 
1.19 
1.19 
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Figure 5. Plot of the binary diffusion coefficient, D\-i, of toluene in supercritical 
carbon dioxide at 313.83 K versus the density of supercritical carbon dioxide. 
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Appendix 1 

Some important hydrodynamic parameters: 



1. 



2. 



Radius ratio: 



Reynolds 
number: 
Re=2ao Uop/t] 



where R^ is the overall 
radius of the coil, and Oq is 
the internal radius of the 
diffusion tube. The effect 
of diffusion tube coiling 
may be considered negli- 
gible if the radius ratio is 
greater than 100. 



where Mq is the average 
carrier fluid velocity, p is 
the carrier density, and rj 
is the carrier viscosity. 
The fluid flow is consid- 
ered laminar if the 
Reynolds number is less 
than 2000. 



3. Schmidt 
number: 
Sc =iq/p Di2 



where D12 is the binary 
diffusion coefficient. 



4. Dean number: 
De=Re Qi-^^ 



The effects of diffusion 
tube coiling can be con- 
sidered negligible if 
ScDe' <20. 
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